Fermi Velocity Spectrum and Incipient Magnetism in TiBe 2 
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We address the origin of the incipient magnetism in TiBe2 through precise first principles calcula- 
tions, which overestimate the ferromagnetic tendency and therefore require correction to account for 
spin fluctuations. TiBe2 has sharp fine structure in its electronic density of states, with a van Hove 
singularity only 3 meV above the Fermi level. Similarly to the isovalent weak ferromagnet ZrZn2, it 
is flat bands along the K-W-U lines of hexagonal face of the fee Brillouin zone make the system prone 
to magnetism, and more so if electrons are added. We find that the Moriya B coefficient (multiply- 
ing ^ in the fluctuation susceptibility Ax(?, w)) is divergent when the velocity vanishes at a point 
on the Fermi surface, which is very close (3 meV) to occurring in TiBe2. In exploring how the FM 
instability (the q=0 Stoner enhancement is S ~ 60) might be suppressed by fluctuations in TiBe2, 
we calculate that the Moriya A coefficient (of q 2 ) is negative, so q=0 is not the primary instability. 
Explicit calculation of Xa(q) shows that its maximum occurs at the X point (1, 0,0)^; TiBe2 is thus 
an incipient antaferromagnet rather than ferromagnet as has been supposed. We further show that 
simple temperature smearing of the peak accounts for most of the temperature dependence of the 
susceptibility, which previously had been attributed to local moments (via a Curie- Weiss fit), and 
that energy dependence of the density of states also strongly affects the magnetic field variation of 
X- 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The cubic Laves compound TiBe2 was already shown 
forty years ago to have quite unusual behavior of the 
magnetic susceptibility x(T) and the Knight shift. 0| x _1 
showed a strong increase with lowering temperature but a 
clear deviation from Curie- Weiss form, while the Knight 
shift was temperature dependent and negative. The mag- 
netic properties of TiBe2 have been controversial since 
Matthias et aZ.0 interpreted the susceptibility peak at 10 
K in TiBe2 as itinerant antiferromagnetism (AFM) with 
an associated magnetic moment of 1.64/is, and Stewart 
et al. reported a transition at 2 K that seemed charac- 
teristic of magnetic ordering. 

However, a clear picture has emerged gradually after 
the idea of weak itinerant antiferromagnetism had been 
abandoned because of the subsequent lack of experimen- 
tal evidence^, Q. Many experiments have shown that 
TiBe2 is instead a strongly enhanced paramagnetic, IE 
and undergoes a metamagnetic transition^ M. Iiofffield- 
driven ferromagnetism) around 5.5 T. Also one can see 
similarity to the magnetic behavior of Ni3Ga by compar- 
ing the values of the low temperature susceptibility, \ = 
1.65 x 1CT 2 emu/mole for Ni 3 Ga[ll| and X = 0.90 x 10" 2 
emu/mole for TiBe2|2j. Based on the magnetization data 
of Monod et aZ|6j Wohlfarth 9j] suggested the transition 
at 5.5 T should be first order. Wohlfarth's considerations 
received at least partial support from theoretical band- 
structure considerations coupled with the de Haas-van 
Alphen data of van Deursen et q?[T3|. 

Clarity began to arise with the extensive experiments 
of Acker et al. who interpreted their magnetization 
datajfj in fields to 21T and the magnetization data of 
Monod et al. @ as evidence for exchange-enhanced para- 



magnetism or spin fluctuations in TiBc2. They found 
the system TiBe2- x Cu x to become FM at a critical con- 
centration x cr — 0.155. Stewart et aL|l3j] measured the 
specific heat of TiBe2 (7 = 42 mJ/mole K 2 ) at low tem- 
perature in and 7T and interpreted the behavior as 
evidence of spin fluctuations. 

The isoelectronic isostructural material ZrZn2 is con- 
sidered a classic example of an weak itinerant ferromag- 
net. Magnetic measurements find very small magnetic 
moments (values from 0.12 to 0.23 \x B )0 El, hence 
the characterization as a weak ferromagnet. The magne- 
tization of ZrZn2 increases substantially with field, but 
unlike TiBe2 with its metamagnetic transition, the in- 
crease continues smoothly to fields as high as 35 T. The 
Curie temperature Tc drops approximately linearly with 
pressure, from 29 K at P = to 4K at P = 16 kbar, 
which extrapolates to a quantum critical point (QCP) at 
P = 18 — 20 kbar. The report of superconductivity coex- 
isting with ferromagnetism in ZrZn2 near this QCP[l6| 
enlivened both theoretical and experimental attention, 
but more recently it has been shown [T^ there is no bulk 
superconductivity. TiBe2, on the other hand, has been 
nearly addressed only rarely for the past twenty years. 

The complex temperature-field behavior of TiBe2 has 
led to many speculations about the microscopic mecha- 
nisms. Of course spin fluctuations play a central part, 
and the highly enhanced susceptibility suggests this sys- 
tem is near a quantum critical point (at slightly en- 
larged lattice constant, say, as well as for the Cu al- 
loying). If FM fluctuations dominate, then a metam- 
agnetic transition (field-driven FM state) around 5 T 
would make sense. If AFM fluctuations dominate, ap- 
plication of a field suppresses the fluctuations, provid- 
ing another way to interpret specific heat under ap- 
plied field. [l8j The anomalies in the conduction elec- 
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tron spin resonance (CESR) linewidth[l^| around 2 K 
have been interpreted in terms of a thermal spontaneous 
magnetism. |20j and a decrease in the resistivity is also 
seen at that temperature. 5] All of these scenarios are sen- 
sitive to the Fermi surface shape, velocity spectrum, and 
possibly the energy dependence of the density of states 
near the Fermi energy, and it is these questions that we 
address in this paper. 

Band structure intricacies by themselves also can come 
into play. Shimizu showed [2(ij that an independent elec- 
tron system with magnetic coupling can undergo a first- 
order transition to a "spontaneous thermal magnetism" 
state (within a range T\ < T < T2) if it is highly en- 
hanced and if the Fermi level lies within a local minimum 
in the density of states. The effects of magnetic fluctu- 
ations should of course be added|2lj to the free energy 
of both the ordered and disordered phases to make this 
treatment more realistic. 

Local density approximation (LDA) energy band stud- 
ies of TiBe2 have been reported previously |U |24| . 
Those studies revealed a split narrow peak in in the 
density of states (DOS) N(E) near the Fermi energy 
(Ep), with calculated Stoner factors IN(Ep) greater 
than unity, giving the Stoner instability to FM. Here / is 
the Stoner exchange interaction averaged over the Fermi 
surface. Thus, as for a few cases that have come to light 
more recently, pH |2(J ferromagnetism is incorrectly pre- 
dicted, indicating the need to account for magnetic fluc- 
tuations not included in LDA that will suppress magnetic 
ordering. By comparing the calculated value of N(Ej?) 
with the measured susceptibility, a Stoner enhancement 
S = [1 - IN(Ep)}- 1 w 60 was obtained, making TiBe 2 a 
more strongly exchange enhanced metal than Pd. 

All of these calculations, carried out 25 years ago, 
used shape approximations for the density and potential, 
and for a detailed investigation of the weak ferromag- 
netism precise electronic structure methods are required. 
In this work, the precise self-consistent full potential 
linearized-augmented-plane-wave (FLAPW) method and 
full potential local orbital minimum basis band structure 
scheme (FPLO) are employed to investigate thoroughly 
the electronic and magnetic properties of TiBe2 based on 
the density functional theory. We compared and checked 
the calculation results of the both methods. We consider 
the effect of magnetism on the band structure and Fermi 
surface, Fermi velocity and compare with experiment and 
previous band calculations. 



II. CRYSTAL STRUCTURE 

TiBe2 crystallizes into a cubic Laves phase C15 crystal 
structure. The C15 (AB 2 ) structure is a close packed 
structure and the site symmetry is high for the two con- 
stituents. Ti atoms occupy the positions of a diamond 
sublattice while the Be atoms form a network of intercon- 
nected tetrahedra, with two formula units per cell. Since 
the major contributions to N(Ep) come from Ti, the lo- 



cal environment of Ti atoms is particularly important to 
keep in mind. Each Ti is surrounded by 12 Be neigh- 
bors at a distance of 2.66 A and tetrahedrally by four Ti 
neighbors a distance 2.78 A away. The TiBe2 structure 
belongs to the Fd3m space group with Ti occupying the 
8a site, and Be the 16d site. The site symmetry of Ti 
is 43m(tetrahedral) and Be has 3m site symmetry. The 
atomic positions are symmetry determined, and we used 
experimental lattice constant 6.426 A for all calculations. 



III. METHOD OF CALCULATIONS 

We have applied the full-potential nonorthogonal local- 
orbital minimum-basis (FPLO) scheme within the local 
density approximation (LDA). [23 In these scalar rela- 
tivistic calculations we used the exchange and correlation 
potential of Perdew and Wang. psj Ti 3s, 3p : 4s, 4p, 3d 
states and Be 2s, 2p, 3d were included as valence states. 
All lower states were treated as core states. We included 
the relatively extended semicore 3s, 3p states of Ti as 
band states because of the considerable overlap of these 
states on nearest neighbors. This overlap would be oth- 
erwise neglected in our FPLO scheme. Be 3d states were 
added to increase the quality of the basis set. The spatial 
extension Of the basis orbitals, controlled by a confining 
potential (r/ro) 4 , was optimized to minimize the total 
energy. 

The self-consistent potentials were carried out on a 
mesh of 50 k points in each direction of the Brillouin 
zone, which corresponds to 3107 k points in the irre- 
ducible zone. A careful sampling of the Brillouin zone 
is necessary to account carefully for the fine structures in 
the density of states near Fermi level Ep. For the more 
delicate numerical integrations, band energies were ex- 
tracted from FPLO in an effective mesh of 360 k points 
in each direction. A separate tool was developed to ex- 
tract energy isosurfaces with gradients from the scaler 
energy grid. The isosurfaces were then used to calculate 
density of states and velocity moments. 

To check carefully the fine structure that we will dis- 
cuss, we also repeated several calculations with the gen- 
eral potential linearized augmented plane wave (LAPW) 
method, as implemented in the WIEN2K code.|3q 
Relativistic effects were included at the scalar relativis- 
tic level. However, we verified that the magnetic mo- 
ment with the experimental structure is not sensitive 
to the inclusion of the spin-orbit interaction. For the 
generalized gradient approximation (GGA) calculations, 
we used the exchange-correlation functional of Perdew, 
Burke, and Ernzerhof. |3J We choose the muffin-tin 
spheres Rmt — 2.6 a.u. for Ti, Rmt — 2.1 a.u. for 
Be and a basis set determined by a plane- wave cutoff of 
RMrK max — 7.0, which gives good convergence. The 
Brillouin zone samplings were done using the special k 
point method with 1280 points in the irreducible zone. 
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FIG. 1: The full LDA band structure of non-magnetic TiBe2 
along symmetry lines showing that there are several bands 
near the Fermi level (taken as the zero of energy) with weak 
dispersion; they are primarily Ti 3d in character. 
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FIG. 2: Band structure of non-magnetic TiBe2 of Fig. 0on 
an expanded scale near Fermi level. The flat bands along L- 
W-U/K-L lines (the hexagonal face of the fee Brillouin zone) 
give rise to the density of states structure discussed in the 
text. 



IV. RESULTS AND DISCUSSIONS 

For orientation we first show the full nonmagnetic band 
structure of TiBe2 in Fig. ^ which is consistent with 
earlier calculations of (23, 123, |24|. The Be 2s bands lie 
between -8 eV and -2 eV. Above them the bands are of 
mixed s,p character, centered on the Be as well as the Ti 
site. Near the Fermi level there are several bands with 
weak dispersion, being of primarily Ti 3d character. The 
bands at K and L are hybridized strongly, while at X 
the s,p character is the main character. As noted also 
by Jarlborg and Freeman. [22l |23^ | one band at L falls ex- 
tremely close to Ef (3 meV below). This band is doubly 
degenerate along T-L, and the L point forms the maxi- 
mum of band 15 and a saddle point for band 16. As the 
Fermi energy rises (for added electrons, say) the Fermi 
surface sweeps through the L point saddle, where the 
band has a vanishing velocity by symmetry. This vanish- 
ing velocity is discussed below. There is another doubly 
degenerate band very near E/ at the W point. 

The density of states (DOS) is shown near Ep in Fig. 
The Fermi energy Ep falls extremely close to the edge 
of a very narrow peak in the DOS. The DOS peak arises 
from Ti d bands hybridized with Be p states. Flat bands 
close to Fermi level centered mostly in regions near the L- 
W-U and W-K directions, i.e. the hexagonal faces of the 
Brillouin zone, cause the sharp peak. Stewart et aL[T^| 
measured the linear specific heat coefficient for TiBe2 of 
7=42 mJ/K 2 mole-formula unit. The calculated value of 
N(Ep)=5.33 states/eV/f.u. for TiBe2 corresponds to a 
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FIG. 3: The total and atom-projected density of states (Ti, 
short dashed line; Be, the lower, long dashed line) for non- 
magnetic TiBe2 per primitive cell. The inset gives the density 
of states for the ferromagnetic TiBe2 showing the exchange 
splitting 0.6 eV. The peak of the DOS for the majority spin is 
entirely below the Fermi level while that of the minority spin 
is above the Fermi level. 
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FIG. 4: (Color online) Fermi surfaces, from top: band 14, X- 
centered pillows; band 15, primarily X-centered jungle gym; 
band 16, F-centered pseudocube; band 17, F-centered sphere. 
Fermi velocities colored dark (red) for lowest to lighter (blue) 



bare value 7 Q =12.6 mJ/K 2 mole(formula unit), leading 
to a thermal mass enhancement 1+A=3.3, or A=2.3 aris- 
ing from phonons, magnetic fluctuations, and Coulomb 
interactions. 

Density functional calculations are usually reliable in 
calculating the instability to ferromagnetism. The en- 
hanced susceptibility|32j is given by 

X{T) = 1-N(E F )I 55 SX0 - (1) 

where xo = ^^(Ep) is the bare susceptibility obtained 
directly from the band structure and / is the Stoner ex- 
change interaction constant. Here N(E^) refers to both 
spins, and hence forward we quote susceptibility in units 
where \xb = 1- The calculation of I is from fixed spin 
moment calculations 33j , in which the energy E(m) is 
calculated subject to the moment being constrained to 
be m. The behavior at small m is E(m) = (l/2)x~ 1 m 2 
from which I = 0.22 eV can be extracted from Eq. ^ 
This value of / gives IN(Ep) = 1.2, larger than unity 
and very close to that calculated earlier, |23| correspond- 
ing to a Stoner ferromagnetic instability. 

As for a few other compounds, TiBe2 is incorrectly pre- 
dicted by LDA to be ferromagnetic. Since spin-orbit cou- 
pling is small in 3d magnets, we neglect it, so the direc- 
tion of magnetic polarization is not coupled to the lattice. 
We have calculated a consistent magnetic moment for 
TiBe 2 : 0.97At B /f.u.(FPLO, LDA), 1.00£t B /f.u.(LAPW, 
LDA), 1.10/z B /f.u.(LAPW, GGA). This value is consid- 
erably larger than an earlier calculation[22j (which also 
reported a much smaller value for ZrZn 2 than obtained 
from more recent calculations!^). We address the over- 
estimate of the tendency to magnetism below. 



A. Fermi Surface and Fermi Velocity 

In Fig. 01 we show the nonmagnetic Fermi surfaces 
shaded by the Fermi velocities. The position of Ep near 
L and W points sensitively determine the exact shape of 
some Fermi surfaces. The shapes can be characterized as 
(a) small F-centered electron sphere from band 17, (b) 
large T-centered electron pseudocube from band 16, (c) 
multiply connected surface mostly enclosing holes around 
the X point from band 15, which we refer to as the jungle 
gym, and (d) flat hole pillows centered at each of the 
three X points. The doubly degenerate bands crossing 
Ep along T-X and X-W guarantee touching of certain 
surfaces along these lines. 

The DOS peak at and above E^ is due to the band near 
the L point where the cube-shaped surfaces are about 
to form bridging necks. Figure shows how the Fermi 
velocity spectrum (N(V; E)) changes with energy at the 
peak just above Ep, at Ep, and at the first minimum 
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FIG. 5: (Color online) Fermi velocity spectrum of TiBe2. The 
low Fermi velocity states are the primary source of changes 
to the density of states. 



below Ep. The Fermi velocity spectrum is defined as 



N(V;E) = Y,KE%-E)5(V % -V) (2) 



C(V;E) Wk X Vfc|Ufc||' 
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with normalization / N{V; E)dV = N(E). Rere£{V;E) 
is the line of intersection of the constant energy Ek = E 
surface with the constant velocity surface \vk \ = V. The 
gradient of the velocity in the denominator makes this 
distribution delicate to calculate accurately. N(E, V) 
was calculated numerically by extracting a triangulated 
energy isosurface from the band structure, then obtain- 
ing a velocity histogram of the states associated with the 
isosurface. 

The spectrum in Fig. [3] shows, at , velocities ex- 
tending down to the very low value of 2xl0 6 cm/s, and 
up to 5x 10 7 cm/s, a variation of a factor of 25. Roughly 
half of the weight lies below 10 7 cm/s. At the van Hove 
singularity at +3 meV, the only noticeable difference is 
additional velocities extending down to zero due to the 
vanishing velocity at L (we have not worried about repro- 
ducing the V — * behavior precisely). At -25 meV, which 
is just below the narrow peak at E^, the strong weight 
in the spectrum appears only at 7xl0 6 cm/s. Note that 
there is very little change in the high velocity spectrum 
over small changes in energy. 
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FIG. 6: (Color online) Top panel: < ^pry > plotted versus 
energy, showing the square root divergence of the inverse mo- 
ment of velocity near the Fermi energy. Unit conversion is: 
1 eV Bohr = 8xl0 6 cm/s. Bottom panel: the graph of the 
second moment of velocity (with constants included to show 
it as the square of the Drude plasma energy) is concave down- 
ward, which gives rise to the negative value of the Moriya A 
parameter. This sign of A is verified by the calculation of 
X(q) a t small q (see text). 







V. ANALYSIS OF VELOCITY DISTRIBUTION 
AND SUSCEPTIBILITY 

A. Renormalization due to Spin Fluctuations 

Following the work of Larson, Mazin, and Singh 35] 
for Pd which builds on Moriya theory, we first attempted 
to identify the relevant band characteristics in order to 
evaluate the spin fluctuation reduction of x m TiBe2. For 
this, one begins with the bare susceptibility in the small 
q and small u> limit, given by 



Xo (q, u) = N(E F )[1 A{^-f + 1 1 < 1 >F -], (3) 

2ir 2 v q 

while the screened susceptibility using the RPA approx- 
imation is given by 



.1 1 



(4) 



The Moriya parameter A = —1.8, expressed in dimcn- 
sionless form here, and mean inverse Fermi velocity < 



1/v >i 



(the second Moriya parameter, discussed 



below) are derived from velocity moments and DOS of 
the band structure, and like the density of states, they 
are greatly influenced by the Fermi surface topology and 
its velocity spectrum. Specifically, changes in topology 
which give rise to points of zero velocity in the band 
structure near the Fermi surface become an important 
factor. The mean inverse Fermi velocity which governs 
the imaginary part of xo(q,w) is given by 



1 



v(E) ^ 



S(s k - E) 



/J2S(e k -E) (5) 



evaluated at Ef- The difference between < v^ 1 >f and 
l/< v >f is one measure of the velocity variation of the 
Fermi surface. The bottom or top of a three-dimensional 
band (corresponding to the appearance or vanishing of 
a Fermi surface) gives only a discontinuity proportional 
to the square of the band mass. At a saddle point, such 
as the merging of the corners of the pseudocube Fermi 
surfaces, v (E) undergoes a 1/ \JE — E cr divergence be- 
cause the associated Fermi surface area does not vanish. 
This "van Hove singularity" in v^ 1 (E) is evident for the 
band edge 3 meV from E^ in TiBc 2 in Fig. |SJ We cal- 
culated 1/vp 1 = 5 x 10 6 cm/s for TiBe 2 . 

For cubic structures, the parameter A in Eq. |3]is given 
by 



A = 



1 



nl(E F ) 
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27T 2 d 2 & p (E F ) 

2 n 



dEp 



(6) 
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N(E F )v 
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Thus A it is proportional to the second derivative of the 
square of the Drude plasma energy il p (i.e. h is ab- 
sorbed into £l p , so flp here explicitly has energy units; 
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FIG. 7: (Color online) Intraband contribution to the real part 
°f x(o)- The increase at small q confirms the sign of Moriya A 
coefficient (see text). Although both [110] and [111] directions 
have a maximum at the zone boundary, the peak along [100] 
(X point of the zone) dominates the instability. 



k sums are understood to be normalized over the zone). 
The second moment of velocity is finite everywhere, but 
its second derivative is not (for example, for free elec- 
trons this diverges as the band edge). Derivatives have 
the unfortunate property of amplifying noise in numerical 
evaluations. We have addressed the noise issue by using 
a large number of k points in the numerical integration 
(360x 360x 360). By fitting VL p (E) 2 with a polynomial 
near the Fermi energy, we obtain the above-mentioned 
value A = —1.8. The Fermi velocity was calculated to be 
v F = 2.3 eV bohr = 1.8 xlO 7 cm/s. 



B. q-dependent Susceptibility 



The negative value of the A parameter indicates, from 
Eq. |3 that the primary magnetic instability in TiBc2 
does not lie at q=0 but rather at finite q, so it is more 
susceptible to AF instability (including possibly a spin 
spiral) rather than ferromagnetic. The sign of A has 
been verified independently by explicit calculation of the 
real part of xio)> with results shown in Fig. \7\ 

The calculation of Xa,p{<i) between bands a and /? 
was performed by an isosurface slicing method. The 
susceptibility can be written, after inserting a factor 
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C. Temperature Dependence of Susceptibility 



A ' 



(7) 



where N*JA) is a susceptibility density calculated from 
the isosurface defined by the Fermi functions and the 
energy S function. The Brillouin zone was divided into 
a 140 x 140 x 140 grid of cubes. Within each cube the 
A integral is calculated as a discrete sum, using variable 
step sizes in A corresponding to 1/30 of the maximum 
difference in energies e„ r -+ e r within the cube. 

The susceptibility rises equally along all three symme- 
try directions (as required by cubic symmetry), but only 
for q along the cubic axis does x{o) continue to increase 
strongly beyond the small-g region. The maximum of 
Xo(q) occurs at the X point, where the intraband part has 
increased by nearly 50% over its q=Q value. In such cases 
where q=0 is not the maximum, it is necessary to apply 
the extension of weak ferromagnets to the AF case.|36j 

The band- by-band contributions to Xo{q) have been 
evaluated both to verify the code and to identify the 
source of the important contributions and structures. 
The sphere FS gives rise to a Lindhard type form with 
2kp « ir/a (but is not perfectly round). The pillows lead 
to a cusp for (q x ,0,0) for q x w 0.287r/a, and along all 
three directions decreases for q > ir/a. For the jungle 
gym and the pseudocube,x increases by a factor of two 
at the zone boundary along (q x , 0,0), with much less vari- 
ation in the other two directions. The contributions to 
N(Ejr) from each of the bands is: sphere, 1.4%; pillows, 
7%; jungle gym, 33%; pseudocube, 58%. 

Away from g=0 the interband contributions to xil) 
contribute, and it is known in other transition metals 
and their compounds that the g-dependence of matrix 
elements can be important. We have calculated also the 
interband xio) f° r several bands around the Fermi level, 
finding that they contribute a broad maximum at inter- 
mediate \q\. It seems unlikely, however, that interband 
contributions will move the maximum away from the X 
point. 

Peaking of xio) a ^ the zone boundary implies a short 
wavelength X — a AF instability (incipient, since no AF 
phase is observed) . With the fee lattice and two Ti atoms 
in the primitive cell, there several possibilities for the 
most unstable mode, which will involve antialignment of 
spins or charge density wave variation, but also may in- 
volve noncollincar alignment of the spins. We have tried 
to obtain a q = AF state within LDA, with atomic 
moments antialigned on the bipartite Ti lattice, but the 
moment vanished when this was tried. We have not in- 
vestigated possible q = X point AF states. 



The high narrow peak in the DOS near E^ suggests 
an explanation of the T-dependence of x mentioned in 
the Introduction, or at least part of it. To understand 
what part arises from simple thermal smearing, we have 
evaluated 



N(E, T) 



df(E - M (r)) 

dE 



}N(E)dE, (8) 



where the chemical potential fJ.(T) is adjusted at each 
temperature to keep the number of electrons (occupied 
states) constant. The result is shown as a series of curves 
for T ranging from zero to 300 K. It is necessary to in- 
clude the variation in fi, and the value of N(fi(T),T) 
decreases by 8%. 

The resulting change in the physical, enhanced suscep- 
tibility is given by 



X(T) 



N(y(T);T) 
l-IN(fx(T);TY 



(9) 



Adjusting I to reproduce the peak height (at 10 K, ex- 
perimentally), which requires 1=0.183 eV (S=56 at the 
maximum of NQj,)), the resulting enhanced Jf(T) is com- 
pared with the data in the lower panel of Fig. |SJ It is 
evident that this simple temperature smearing accounts 
for much of the observed temperature dependence. Ad- 
ditional indirect temperature smearing will come from 
phonons and from electronic and magnetic interactions 
as these excitations are increasingly excited upon raising 
the temperature. We conclude that TiBe2 contains no 
appreciable contribution to the susceptibility from local 
moments. 



D. Field Dependence of Susceptibility 

For an energy-dependent DOS and a highly enhanced 
susceptibility, a field-dependent susceptibility x(H,T — 
0) = x{H) is expected. In TiBe 2 a strong effect of this 
kind has been seen, which can be characterized as field- 
driven ferromagnetism. The differential susceptibility 
Xd{H) = dM{H) / dH where M is given by the difference 
in electron occupations n a (H). A many-body treatment 
shows that the spin imbalance can be expressed |38fl in 
terms of the spin-dependent thermal (energy E surface 
averaged) Green's function 



G„(E,iu n ;H) 



1 



iuj n - (E - [i- afi B H) - T, a (H) 



(10) 



i AH) = J dEN(E)Tj2GAE,iuj n ;H) 



here w„ is the fermionic Matsubara discrete energy vari- 
able and rj is a positive infinitesimal. The simplest form 
of (Stoner) self-energy = aK^sH should be appro- 
priate (l + K = S). 
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FIG. 8: (Color online) The upper graph shows how the density 
of states near the Fermi energy changes from T=0 to T=300K. 
The lower graph shows the experimental susceptibility |37| 
compared to theory. The Stoner I has been adjusted slightly 
from the calculated value to match the susceptibility maxi- 
mum. 
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FIG. 9: (Color online) Magnetic field (H) dependence of the 
Fermi level density of states for TiBe2, referred to its H=0 
value, as defined in the text. The initial increase with field 
indicates an increasing instability towards ferromagnetic or- 
der. 

dependence of S due to the structure in N(E) and the 
delicate situation here that IN(E F ) is approaching unity, 
but at this point we neglect such details. 
The result for the relative correction 

_ N(E F - SijlbH) + N{E F + Sfi B H) 
K ~ 2N(E F ) 

is shown in Fig. 03 The effect on the ratio (thus on the 
differential susceptibility) is clear, however even with the 
factor of S=60 enhancement of the energy scale {(i B H — > 
Sfi B H) the peak occurs at a field one order of magnitude 
smaller than seen in experiment. This difference seems to 
indicate that the field influence on the spin fluctuations 
dominates; however the variation in N(E) will need to be 
accounted for in any quantitative explanation. 

VI. SUMMARY 



Taking the field derivative of M(H) but keeping H 
finite, and using (in this approximation) 

G a (E, iu n ; H) = G (E- a^ B (l + K)H,iu n ;H = Ojll) 

we obtain the result at zero temperature 

dM{H) 



Xd{H) 



d^ B H) 

= S[N(E F - Sn B H) + N(E F + Sfi B H)] 



(12) 



This clearly reduces to the usual T=0 result at H=0. 
A slightly better treatment would have also some H- 



The complex and sometimes confusing data on the en- 
hanced paramagnet were discussed in the Introduction. 
It seems clear that magnetic fluctuations will be required 
to understand the underlying mechanisms. Here we have 
presented a precise calculation and analysis of the elec- 
tronic structure, especially focusing on the Fermi sur- 
faces and velocity spectrum at and near the Fermi level 
that underlies not only the single particle excitations but 
also the spectrum of magnetic fluctuation in the itinerant 
limit, which clearly seems to be the case in TiBe2. 

Our calculations have confirmed the sharp structure in 
the density of states around the Fermi level that had been 
noted earlier, and quantified the tiny energy scale that is 
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involved: the Fermi level lies in a region of steep DOS, 
just 3 meV from an abrupt van Hove singularity. This 
singularity is derived from a doubly degenerate band at 
the L point of the zone. We have shown how to calcu- 
lated the spectrum of velocities (speeds) over the Fermi 
surface, and find the spectrum to be peaked at (the low 
value of) 10 7 cm/s, with much of the weight below that 
value. Moriya theory for weak ferromagnets requires, for 
the imaginary part of the inverse susceptibility, the mo- 
ment < l/vp >; we have illustrated that X/v{E) diverges 
at the van Hove singularity signaling possible problems 
with applying Moriya theory to TiBe2- 

Moriya theory for weak ferromagnets also requires the 
dimensionless quantity A oc d 2 il p (E)/dE 2 at the Fermi 
energy, where Q p is the conventional transport Drude en- 
ergy. We find that this quantity is not positive, as it must 
be for an incipient ferromagnet; rather it is negative in- 
dicating the dominating (nearby) magnetic instability is 
finite q: antifcrromagnetic, spin wave, spin spiral, etc. 
Direct calculation of the generalized susceptibility Xo(q) 
confirms the sign of A, and reveals the dominant insta- 
bility to lie at the X point of the Brillouin zone, making 
TiBe2 an incipient antiferromagnet. 

We have shown that the sharp structure in N(E) has 
other consequences. First, it leads to a T-dependent 
chemical potential. Together with the temperature 
broadening of N(E) and the Stoner enhancement S « 60, 
this simple temperature broadening can account for most 
if not all of the temperature dependence of the suscepti- 
bility, which some investigators had interpreted as Curie- 



Weiss-like. As a result, the occurrence of local moments 
in TiBe2 can be ruled out. Similarly, we have shown 
that this sharp structure in N(E), again together with 
the large Stoner enhancement, has a substantial effect 
on the field-dependence of the differential susceptibility. 
There is still the question of how much of the measured 
field dependence is due to this induced exchange split- 
ting, and how much is due to the effect of the field on the 
magnetic fluctuations. 

Many of the results we have obtained here are strongly 
dependent on details of the band structure and the po- 
sition of the Fermi level. That these results reflect real- 
istically the mechanisms underlying the many fascinat- 
ing observations obviously requires that the band struc- 
ture formalism is applicable in detail to such systems and 
that the calculations are accurate. Another requirement 
is that of high sample quality, that the stoichiometry is 
precise and that defect concentration must be very low 
(simple impurity broadening will affect behavior). These 
questions must be addressed in deciding whether to press 
onward to a more complete and more challenging expla- 
nation that includes effects of both magnetic fluctuations 
and the energy dependence of the density of states. 
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